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SUMMARY 


A nonlinear simulation model of the CH-47B helicopter, developed by the Boeing 
- 1 J om P an y ( r ef. 1). has been adapted for use in the NASA Ames Research Center 
(ARC) simulation facility. The model represents the specific configuration of the ARC 
variable stability CH-47B helicopter (fig. 1) and will be used in ground ablation 
research and to expedite and verify flight experiment design. 

Modeling of the helicopter uses a total force approach in six rigid body degrees 
of freedom. Rotor dynamics are simulated using the Wheatley-Bailey equations, 
including steady-state flapping dynamics. Also included in the model is the option 
tor simulation or external suspension, slung-load equations of motion. 

n1 . Validation of the model (discussed in Volume II of this report) has been accom- 
moHpl 6 j S f , ? atatic and dynamic data from the original Boeing Vertol mathematical 
model I * 1 aght test data from references 2 and 3, as reproduced in reference 4. The 

ARC SiiL a ?? r ° P l USeiU real ’ tiine P iloted a i">olation and is implemented on the 

ARC Sigma IX computer where it may be operated with a digital cycle time of 0.03 sec. 


NOMENCLATURE 


AERO 

ARC 

BV 

c.g. 

CONTROL 

DCPT 

ECS 

ENGINE 

n b 

rpm 

ROTOR 

SAS 

SLING 


fuselage aerodynamics subroutine 
Ames Research Center 
Boeing Vertol Company 
center of gravity 

mechanical control system subroutine 
differential collective pitch trim 
electronic control system 
engine and governor subroutine 

change in helicopter yawing moment per sideslip angle 

revolutions per minute 

rotor dynamics subroutine 

stability augmentation system 

sling load dynamics subroutine 


SNP 


shaft-normal -plane 


SNPW shaft-normal -plane-wind 

V equivalent velocity 


INTRODUCTION 


At Ames Research Center (ARC) , the CH-47B provides a unique capability for generic 
flight research in flight controls and displays for rotorcraft and VTOL aircraft. In 
addition to the existing potential for variable-stability flight, a programmable dis- 
play system and a variable force-feel system are being developed. The purpose of this 
mathematical model development is to provide the capability for real-time simulation 
and for the preliminary check-out of in-flight research experiments for the variable- 
stability CH-47B helicopter. 

Subroutines that comprise the mathematical model describe the rotor systems, 
fuselage aerodynamics, engine and governor, mechanical control system, the option for 
either an electronic control system or the basic stability augmentation system (SAS), 
and the option for externally suspended, slung-load dynamics. Forward and rear rotor 
dynamics are simulated in a shaft-normal-plane-wind (SNPW) reference frame with the 
Wheatley-Bailey (modified tip path plane) equations of references 5, 6, and 7. Steady 
state flapping dynamics are represented with these equations; however, in-plane 
motions are neglected. Forces and moments at the rotor hubs are then calculated as a 
function of rotor aerodynamic conditions and dynamics, after which they are resolved 
to the helicopter center of gravity. Six rigid-body forces and moments resulting 
from fuselage aerodynamics are found from tabular data interpolated as a function of 
fuselage angle of attack and sideslip angle. 

Each engine is represented with nonlinear, second-order dynamics; left and right 
engine models are identical, yet are modeled separately. The fuel control system and 
gas generator are each modeled as a first-order system, the latter including a variable 
time constant dependent upon power and power error. The engine governor, whose pur- 
pose is to regulate rotor rpm, is modeled as a linear, third-order system. 

Modeling of the hardware from the cockpit controls to the swashplate comprises 
the mechanical controls subroutine. Included are upstream limiters on each control 
input, first- and second-stage mixing, swashplate limits, and swiveling and pivoting 
actuation dynamics (first order). 

Stability augmentation in the form of longitudinal, lateral, and directional rate 
damping is modeled. Additional features of the directional SAS include turn coordina- 
tion and feedback of sideslip angle to obtain a stable yawing moment change with 
sideslip (Ng). 

The provision for an electronic control system (ECS) model has been included in 
this program. Although no specific ECS configuration has been documented in this 
report, the information necessary to integrate such a subroutine into the simulation 
model is discussed in the section concerning the ECS. 

A model of an externally suspended, slung load has been developed and is available 
for use with the helicopter simulation model. Three state variables, defining the 
position ot the load and suspension cables relative to the helicopter, are represented 
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copter and slung^d) , "s^epKsented |" < '“ 0n ' Thus • the combined system (hel i- 
modeliug the two rigid bodies " lne COUpl<! ‘ 1 ' differential equations 

-r are *> 

by an engineering explanation, input/output variabT Vu Subroutlne is characterized 
computer mnemonics in terms of engineering vlr^M 8 f S ’ and the def inition of 
n the following order: rotor dynamics (ROTOR) f^'i ^ subroutines are discussed 

and governor (ENGINE), mechanical control system (COOTRn^ aerodynainics (AERO) , engine 

system (SAS), electronic control system (ECS) and ’ Stabilit y augmentation 

system (ECS) and slung-load dynamics (SLING). 

Operational considerations are discussed dn.i ,1 
constants and other information necessary f 'a ^ ^specification of input 

cab and a visual display, 7 Piloted simulation using a simulator 

“ «P S| ; ««“ »d dynamic mode! 

also ARC dynamic data are compared £ lt h a Boeing' SaS " eriVative da “ are tabulated 
Ight-test data from references 2 and 3 (reproduced in ref” P 4)! <BV) "°‘ lel and CH - 47 

MATHEMATICAL MODEL DESCRIPTION 


Rotors 

basis for'the si^lLio^o^t^^otor^'n 1 this eq “? tlons (refs ‘ 5 ’ 6 * ^nd 7) form the 
ROTOR, total forces and moments resulting fJom m ° del - In su broutine 

the SNPW reference frame. These are then * acb helicopter rotor are computed in 
the helicopter center of gravity (c.e.) for^ 115 ° rmed to the bod y reference frame at 
freedom rigid-body equations (which are part ^ 0 ^ orat10 ^ into the six degree-of- 
facility and are known as subroutine SMART (refs. 8 and^))?^ ^ simulati °n 

putcr" e :L e „L s innc\uiKrv«i o rbu a L r ™:Vd he rotor subrouti - <*■ .« =o»- 

routines. The equations are executed sequential] 0 ^ 1 ^ a° ^ frOTn ° ther 1110(161 sub ~ 

modules in the figure. Since the calc.l^M * 38 indlcated by the numbered 

required for this model, it is necessa * tlon ° f rotor hub forces and moments is 

copter body and the SNPW reference frames *>«“«" the heli- 

totor c g. relative to the actua he lender t° 5 3 ' ^ Portion of the actual 
tlon (O. encopter c.g. (fig. 3 ) ls computed us in g ~ ^ ua - 
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c.g. are given by: ~ ° r C,g * re lative to the baseline helicopt 


er 


3 


% * 


The vector 


20.43 ft 


0.0 ft L J d 


7.49 ft 


DXCG (AX 


DYCG J AY 


-18.46 ft 
0.0 ft 
12. 16 ft 


DZCG AZ 


is the position in inches of the c. g. of the octuai helicopter relative to the base- 
line specifications. Baseline helicopter c.g. positions are 

r x (331 in. 

C H 

J Y - J 0.0 in. > 

I C ’ 8 ' 

Z 11.2 in. 

L c -gJ L J 

and the sign conventions are as given in figure 4. 

To compute forces and moments at the rotor hub, helicopter body-axis velocities 
(from subroutine SMART, rigid-body dynamics model) are transformed from the body 
reference frame to the rotor SNPW reference frame. Representation of the body axis 
velocities at the rotor hubs is given in equation (2). 


UFR1 ,URRl u t ,u 0 
F 1 Ki 


VFR1 , VRRl v ,v D = v + h 

F i R i S 


WFR1 ,WRRl w ,w D 

F i R i 


0 

_h F,R 

d F,R 

f,r 

0 

0 

F,R 

F,R 

-0 

F,R 
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Body-axis velocities (at the rotor hub) are transformed (eq. (3)) from the body to 
the SNP reference frame through shaft incidence angles ip,i£ (fig* 5). 


UFR2 ,URR2 u-, # u D 
F 2 R ? 


VFR2 , VRR2 v r ,v n 
F 2 R 2 


WR2 , WRR2 Iw , w 

I F 2 R 2 


cos i F>R 0 sin i F>R u f ,u Ri 


VX 


-sin i F>R 0 cos i F)R j w F]! w Ri 


The rotor SNP may be considered an intermediate reference frame between the helicopter 
body and SNPW reference frames. 
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Rotor sideslip angle is defined by equation (4), 


BETAFR 


B' = arctan (i 

BETARR ’ F,R 2 

and SNP translational velocities are effectively resolved (eqs. (5) and (6)) through 
Bp R into the SNPW reference frame, as shown in figure 6. Rotor-hub forces and 
moments are eventually computed in this frame, as indicated in the figure. 


U F,R 2 + v f,r 2 


WRR - “F,R 2 

Next, helicopter-body angular velocities (from SMART) are transformed (eqs. (7) 
and (8)) to the SNPW reference frame as shown in figure 6. 


PFR p. 


cos Bp cos ip sin Bp cos Bp sin i p R 


QFR q F = -sin Bp cos i p cos Bp sin Bp sin i p q g 


RFR r. 


PRR p t 


-cos Bp cos ip -sin Bp 


-cos Bp sxn ip P B 


QRR q t 


-sin Bp cos i R cos Bp -sin Bp sin i R q fi 


RRR n 


Rotor angular velocity is corrected for helicopter yaw rate in equation (9) 


OMEGFR 


OMEGRR 


S\R fi F,R r F,R 


and rotor tip speed is calculated based on this rpm in equation (10) 


VTIPFR 


VTIPRR Tlp F,R ^F.R F,R 


Advance ratio and the free stream component of inflow ratio are calculated in equa- 
tions (11) and ( 12) : 
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AMUFR 
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F,R 
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ALMPFR 
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w. 


F,R 
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r F,R^ 


(ID 

( 12 ) 


Prior to their usage in computations (Le. , for flapping coefficients and rotor 
forces and moments), the pilot* s control inputs are transformed to the SNPW reference 
frame and corrected for control phasing angle ( 4>p) and pitch-flap coupling (5 3 ). Thus, 
it is unnecessary to make these corrections during the actual computation of these 
quantities (as noted in the flapping assumptions which follow). Longitudinal and 
lateral cyclic pitch in the SNP reference frame (from subroutine CONTROL) are trans- 
formed to the SNPW reference frame in equations (13) and (14). 
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(13) 


(14) 


Although the pitch-flap coupling and control phasing angles are zero in the current 
configuration of the ARC CH-47B, the capability for these variations has been included 
in the simulation model. The purpose of the control phasing angle, <j>p , is to offset 
the lead of the blade relative to the pitch hinge, which was introduced by pitch-flap 
(6 3 ) coupling (fig. 7, taken from ref. 10). In equations (15) and (16), rotor cyclic 
pitch positions are transformed through control phasing angle, 4?p (fig. 8). 
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(13) 


(16) 


In equation (17), rotor cyclic and collective positions are corrected for 6 3 
(ref. 11). 
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AICFR.AICRR 

A ir 

l f,r 


A !c 

f,r 2 


ai F,R 

BICFR.BICRR 

B l C 

F,R 

as 

B I C 

f,r 2 

+ K ft 
e F,R 

bl F,R 

THOFR.THORR 
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1 

Pn 

— o 
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a °F,R 


( 17 ) 


where Kg F ^ R = -tan(6 3F ^ R ). 

Rotor degrees of freedom are limited to feathering and the computation of steady 
state flapping and coning coefficients. No in-plane (lead-lag) degree of freedom has 
been considered. Flapping and coning coefficients are computed by solving a 3x3 linear 
system of algebraic equations, and are developed based upon the following simplifying 
assumptions (ref. 1): 


1. Only the first harmonic terms are used. 

2. There is a uniform inflow. 

3. No reverse flow is considered. 

4. Identical forms for the front and rear rotors are used. 

5. There are no pitch-flap coupling effects (the control inputs are corrected 
in this regard) . 

6. There is a zero tip-loss factor. 

7. There is a negligible hinge offset. 

8. Rigid blades are used. 

9. There are no compressibility effects. 

10. - There is a constant rotor airfoil-section lift-curve slope. 

11. The rotor airfoil-section drag varies only with rotor angle of attack. 

Steady-state flapping and coning angles are found by solving equation (18) with 
Cramer’s Rule. (The derivation of these equations is given in appendix A.) 
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F.R 4 8 


F F ,R = K e. 
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original r:r ■; 

° F POOR j 


F.R 


f.r 


(*%) 
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’f.r 
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4 - ->q. 
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F.R F * R 


tw (llZ + 1*5 r) 
F.R *’ R 
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a _ Vr\ 

,R C F, rR B Vr \ 4 / 


l f,r ~ A i 


f.r 

161 


C F,R 2 pa 5 


L A* A _ P F. R \ 
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Using Wheat ley- Bailey theory (rpfq s — t\ .1 
the rotor hubs are computed. Expression* f ’ ^ hrust » torque, side force, and drag at 
developed for a tandem^ heUcopter uslnv ^ J° rque follow the theory « ' 
for rotor side force and drag were greatly s I Vlff C ® fraIne • Ex P re ssions 

because the simplified forms providfd a y siraplifi ® d b y Bv during their development 
the full theoretical expre^ions “ tch With fll8ht test data than did 
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Mean rotor thrust is computed with equation (19). 


origins. { ... 
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CTFR1 

CTRR1 




( 19 ) 


In coefficient form, thrust is modified owing to limits on its maximum allowable value, 
for rotor stall, and due to ground effect. Since the maximum allowable normalized 
thrust coefficient, 2Cx/aa, is 1.0, a limit is imposed if the computed value is 
greater than 1.0. As a function of advanced ratio, normalized thrust coefficient is 
modified as shown in figure 9 for the effects of rotor stall. This is an empirical 
correction which was derived by BV to provide a better match of the model's dynamic 
response with wind-tunnel test data and is selected (along with a correction to rotor 
torque) with flag NSTALL in the simulation model. Thrust coefficient is computed as 
shown in equation (20) 


CTFR 


CTFR 


2C 


T F,R \ a F,R°F,R 


F,R 


a F ,R a F,l 


( 20 ) 


and if longitudinal velocity is less than 40 knots (and if the ground-effect correc- 
tion is selected with flag NGREFF) , thrust is modified for ground effect as a function 
of altitude and airspeed. Thrust is calculated in equation (21) 


T = C pirR ' 

F ’ R t f,r b f,r 


2 F r 1 + K e e T i g e 
F,R\ g.e. F , R i.g.e. F 


( 21 ) 


where 
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ge 


= 1 - 


F,R 


F,R 

U 

ge 


(U =40 knots) 

g.e. 


and T. is determined from figure 10 as a function of the rotor height to diameter 

i.g.e. 

ratio (h/D) 

v ' rotor 
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Mean aerodynamic: torque required is found from equation (22) 


CQFRl 2C Q F-R 
CQRR] a F,R°F,R 



'l 

3a i 


F,R 
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1 C F,R blF * R 




+ 






C F,R lp,R 


- a 


1 F,R 



As a function of rotor thrust and advance ratio, the torque coefficient is modified 
for rotor sr.all (flag NSTALL) as shown in figure 11. A] so, an empirical correction 
is made to the torque coefficient to attain a better match with flight-test data. 

This correction, the effects of which are shown in figure 12, is calculated as a func- 
tion of advance ratio and thrust coefficient (flag NTRQCR) . Including the two 
corrections, the aerodynamic torque coefficient is: 


CQFR 

C Q 

CQRR V F,R 


. (*%. rVf/f. 

\ a F , R°F , R J 1 


R 


+ ac q f , + ac q 

* >r r.s. 


(23) 


F,R 


and the rotor torque required is 

QAERFR 

QAERRR 


aer f,i 


C_ pttR^ !>'2 2 

Q F ,r b f,r 


(24 


Rotor sideforce is calculated with equates (25) - (27) 
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2C n 
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,R ^4 ( bl F,R " Al C F>R ) ^F ,R a °F ,R^ 


+ 2 a °F 


5 °°F,R 3X F > R ■ °tw 


+ 4 X F 


- l - 5 

•,r( |51 F,R ■ Ai C f>k ) + 6 a °F,R( Bl C KjR + ai F,R^ 


(25) 


CYFR 
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Y F,R ] a F,R°F,R 
CYRR Y F,R \ a F,R CT F,R/ 2 


(26) 
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YFR 


Y p r " c v P*Ru U v d 
YRR ,R Y F,R “f,r F * R 


( 27 ) 


drac IH^ceW !ffT J® f 8 ® 01 " 641 / 0 ^ blade-profile drag (eq. (28)) and the normalized- 
arag (H force) coefficient is calculated as xn equation (29). 


DELFR 

DELRR 6F ' R ‘ 40r . R + 9Sl F.R C T 


F,R 


(28) 


CHFR1 2C H 
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= C„ 


CHRR1 a F,R"F,R T F,R 
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(29) 


respectively^ ^ <M) Sh °” ^ CalCulatlons for ™tor d ™8 coefficient and !i-force, 
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CHRR H F,R \ a F,R°F ,Ry 
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HFR 


R F R P^RtC d 

HRR F,R H F,R B f,R F ’ R 


(31) 


are found 


ouid^f Ti Pl ^ hin8 r m0inen ^ S at the rotor hub resulting from aerodynamic forces 
ound as a function of steady-state flapping angles (eqs . (32) and (33)). 
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(33) 


Inflow ratio dynamics, which are modeled using the ARC 1 -f ri o a r* -f * 

gram LOLIN frpf }'?\ t g Lne aw- local linearization pro- 
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f \ 2 Y4 tt ! 2 /^ i ; 


(34) 




R fi„R. 


D F C T 
FR l ? 
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Referring to equations (34) and (35), rotor-on-rotor interference parameters 

(rear on forward) and Dr™ (forward on rear) are calculated as shown in equation (36) 
v r R 

BDFFR 


BDFRF 


D ■ d p ( l l sin e FUS^ + C F 2 l Sin ^FUS 

rR(RF) FR(RF) 2 


(36) 


where and dp Rp are found, depending upon whether the helicopter is in forward 

or rearward flight, in figures 13 and 14 and Cp^ is found in 1 igure 15. A more 

detailed description of the LOLIN approach to solving equations (34) and (35) and an 
explanation of the differences between this approach and the one used originally by 
BV, is given in appendix B. 

Finally, rotor forces and moments at each rotor hub are transformed to the heli- 
copter c.g. These forces and moments form a portion of the total forces and moments 

acting on the rigid body (helicopter) and are integrated in SMART to give the trans- 

lational. and rotational states. 

In equations (37) and (38), forces at each rotor hub are transformed from the 
SNPW to the helicopter body reference frame. 
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-cos g R 


sin 8 R sin i R -cos ij. 
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(38) 


Total moments at the helicopter c.g. due to the rotors have contributions from 
two sources: (1) moments at the helicopter c.g. resulting from forces at the rotor 

hub and (2) moments at the hub transformed to the c.g. Equation (39) shows the com- 
putation of the first contribution, equations (40) and (41) show the computation of 
the second contribution, and equation (42) gives the summation of moments from each of 
the two sources. 
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Total rotor forces and moments, 

r XAERFR ,XAERRR^ fALARFR VLARRr’] 

< YAERFR , YAERRR > and < AMARFR , AMARRR V 

1.zaerfr,zaerrrJ [anarfr.anarrr J 


( 40 ) 


(41) 


(42) 


are passed to the AERO subroutine for summation with the fuselage quantities calcu- 
lated therein; aerodynamic forces and moments (rotor + fuselage) are transferred to 
SMART as inputs ' 



Table 1 is a list of the ROTOR subroutine variables together with constants and 
conversion factors. Included is each variable, its FORTRAN mnemonic, units, common 
location, if applicable, and physical description. Table 2 is a list of the variables 
transferred between ROTOR and other subroutines. 


Fuselage Aerodynamics 

Tabular data from rotor-off wind-tunnel tests provides the basis for fuselage 
aerodynamic forces and moments. These are represented in the helicopter body refer- 
ence frame and are normalized by fuselage dynamic pressure. The data are obtained 
from the function tables by linear interpolation on fuselage angle of attack and 
sideslip angle (figs. 16-21). 
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UKlt, • • i 

of pod;? v . 

To calculate fuselage angle of attack, rotor downwash velocity at the fuselage is 
computed with an empirical expression, and is used to modify vertical velocity 
(eq. (43)1 


(^p ^p)^p^2 

™ PR W B ’ ”B 1+- D- ~ <«) 

F RF 

Using the vertical velocity at the fuselage, Wg, fuselage angle of attack is calcu- 
lated from equation (44). 


ALPHFS 



(44) 


Fuselage sideslip angle is computed in equation (45) , which is somewhat simplified 
from the helicopter sideslip angle computed in SMART. 


BETAFS 


Bpug = arctan 



(45) 


Fuselage dynamic pressure, used to normalize force and moment entries in the function 
tables, is found using equation (46). 


SQFS 


’FUS 


P(u 2 + 


+ < 2 ) 


(46) 


From the function tables, the resulting forces and moments are: (D/q)pyg, 

( Y / q ) FU S > (L/q)pug, CSf/q) p-yg , (M/q)pyg, (N/q)pyg. These quantities are then corrected 
for differences in the equivalent "flat-plate area" between the actual helicopter and 
the model used in the wind-tunnel tests from which the data were obtained. This 
correction accounts for additional sources of drag (i.e., rotor hubs, rotor blades, 
landing gear) that were not included in the wind-tunnel model. 


Correction terms to be applied to the fuselage forces are calculated as shown in 
equations (47)-(49), where Afe is the difference in flat-plate area; fuselage forces 
are calculated in equations (50)-(52). 



Afe 


[1 + 

(tan a pus ) (tan $ FUS ) 2 

] 1/2 


Afe tan Bp us 


[1 + 

(tan a FUS ) 2 (tan B Fug ) 2 ' 

I l/2 


Afe tan 



[! + (tan a FUS ) 2 (tan 


XAERFS 



(47) 


(48) 


(49) 


(50) 
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YAERFS y fus = q FUS 


ZAERFS Zpyg - q pus 


[0L - •©„,] 
fcL ‘ 


(51) 


(52) 


To make the corrections necessary for differences in c.g. position between the actual 
helicopter and the wind-tunnel model, this moment arm is computed as in equation (53) 


SLCFS 

V. 

C FUS 


h 1 

c 

X 


"ax / 1 2 

c.g. 


SDCFS 

d 

C FUS 

= 

d 

c 

X 

- 

AY / 12 

c.g. 

(53) 

SHCFS 

h 

. C FUS_ 


o 

1 


AZ / 1 2 

c.g. 



\ ' 

x 


-1.47 

ft' 


/X 


is the baseline model c.g. offset (fig. 22), which 




c.g. 

d^ 

0 ft 



/Y 

c.g. 

c x 

has the constant numerical value of 


• 

In 


1.31 

ft 


Z 

L c >J 



m 


c.g : 


is 


«■ 

li- 


the position (in inches) of the c.g. of the actual helicopter relative to its baseline 
(fig. 4). 

Using equation (54), fuselage moments are adjusted for this difference in c.g. 
position , 


ALARFS 

AMARFS 

ANARFS 


l fus 


( y, / q ) F us 


^us 

= 

( : */q) fus 

q FUS + 

n fus 


( N/q ) F us 



0 


-h 


'F US 
-d 

L C FUS 


'FUS 

0 

'FUS 


-a 


'FUS 

'FUS 

0 


FUS 


FUS 


'FUS 


(54) 


If the helicopter is in rearward flight, the signs on Xpus* M FUS» and N FUS are 
reversed to account for the aerodynamic differences at this flight condition. 

Total aerodynamic forces and moments include rotor and fuselage contributions, 
which are summed at the end of the subroutine (eqs. (55) and (56)) and passed to 
SMART. 


FAX 

X AERO 


X FUS 


pt 

w 

< 




Xaer r 

FAY 

y aero 

- 

y fus 

+ 

y aer f 

+ 

yaer r 

FAZ 

Z AER() 

m m 


Z FUS 

» m 


z aer f 

■» m 




( 55 ) 
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L 

AERO 


FUS 

TAM 

m aero 

= 

^US 

TAN 

n aero 


n fus 


L AER f 


1 ' aer r' 

m aer f 

+ 

m aer r 

n aer f 


***** 


( 56 ) 


Table 3 gives the definition of the variables, constants, 
of the AERO subroutine. Table 4 lists input /output variables 
routines, together with required input data. 


and conversion factors 
to and from other sub- 


Engine and Governor 

P ° wer Js supplied to the rotor system by two Lycoming T-55-L7C turbine engines 

each ensinP 1* J 77*' Alth ° ugh the representations are identical mathematically, 
engine is modeled separately. The block diagram in figure 23 illustrates t..e 

dvnam! 1 " 8 me J ho< ? for tbe left en 8 ine ’ including the governor and forward rotor-shaft 
dynamics. Nonlinear functions are shown in more detail in figures 24-30. 

. . , Aa sho ™ n ln figure 23, trimming of the engine by zeroing Q , is done while in 
initial condition (I.C.) mode by setting flag ISTEADY after the rigid-body states have 

oTl^rT,' \- Pil0t .T^ S ’ Sh °"" ° n the left Slde ° f tl “ diagram! include poHt^! 
droon^ m 1 ^ atlck (6 - 13 fed forward into the engine to compensate for rpm 

° '’/i l6Ver (compressor speed); and beep trim switches (torque may be adjusted 

e t engine individually, or engine torque and rotor rpm may be adjusted on 

slI " ultaneousl y)- Changes in beep trimmer and collective positions modify 
the fuel control actuator (N 2 ) command. The fuel control mechanism is modeled as a 
l.st order system, with friction in the response represented by a deadband and by 
hysteresis . Unlimited commanded power is calculated, as shown in figure 31, as the 
difference between equivalent rotor rpm (N R ) and the fuel control actuator position 
%; ter "* % provides the intercept of the unlimited commanded power curve, 

and the slope of the^curve (M) is an empirically derived constant between engine f^i 
low and engine power. Feedback of Nr ( ft) in the unlimited, commanded power calcula- 

tion represents the governing loop of the engine, where engine power is modified to 
regulate variations in rotor rpm. ouxiieu to 

tion nf S n°7\ in 7 & 7 eS 23 and 31 ' the to PP ing P° wer level of the engine is a func- 
are available o^the * CO,n P reS8 ° r speed ( N i>- Three positions, STOP, GROUND, and FLY 
stant rate of 0 8 In / 1 actaator motion between the positions is at a con- 

' 8 in -/ sec - Unlimited commanded power is then topped as a function of 
rotor rpm and compressor speed. tunccion ot 

internal 71117^17^77 1 m ° deled 38 3 flrSt -° rder lag with a time constant and 
stare is i r 3r6 variable - The gas generator dynamics time con- 

star.r is a function of power output, modified as a function of power error The 

variable internal limiter adjusts for the engine, which powers Sown six times 
as er t tan it powers up, and is a function of power output. 

ffl. T ^ en *:r, fi r r T an ' 1 rotor shaft dynamic., modeled a, a third-order system 

power available ™ ! "h rP "; to the *° vernor '’ nd shaft dynamics model are 

power available from each engine and power required for the accessories (hydraulic 

hfd'me-rlT 0 " l 0 S T’ As — *" figure, this s^ ei Is dr ten by 

the difference between resistive torque (damping plus spring torque) and rotor torque 
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required. Rotor acceleration is the difference between shaft resistive torque and 
engine torque available. Engine outputs: rot^r rnm (OMEGA), spring torque 

f qgovf! 

|qgovrj 

and rotor rpm uncorrected for helicopter yaw rate 

f omegpfI 

^OMEGPRJ 


are passed to the ROTOR subroutine. 

Table 5 is a list of variables computed in the ENGINE subroutine, together with 
constants and conversion factors, and table 6 has ENGINE subroutine input/output 
variables, and logical flags. 


Mechanical Controls 

The purpose of the CONTROL subroutine is to represent the mechanical hardware 
between the cockpit controls and the rotor swashplate. A block diagram of the sub- 
routine logic is shown in figure 32. Mechanical control system inputs are lateral 
cyclic ( 6 Ap ), collective (5c p ), longitudinal stick (6 b p )> and directional ( 6 Rp ) cock- 
pit control positions. The SAS and ECS actuator inputs from the respective subrou- 
tines, and selected with the flags shown in the figure, augment the appropriate 
cockpit control positions. Longitudinal cyclic position is also augmented by the 
differential-collective-pitch-trim (DCPT) actuator which (although this capability 
has been disconnected in the ARC helicopter) may be selected in the simulation model 
by setting flag IDCPT. The purpose of the DCPT actuator is to artificially provide a 
stable longitudinal stick position gradient with airspeed (fig. 33). To accomplish 
this, as a function of airspeed, the DCPT actuator automatically introduces a positive 
pitching moment (fig. 34), requiring the pilot to move the longitudinal stick forward 
to maintain trim (ref. 13). 

After control-stop limiting (downstream of cockpit control-position limiting, 
which is not included in the diagram) , control positions are converted from inches to 
degrees of equivalent swashplate, resulting in © Ap R , ©c p R » ©Bp R > an d ©Rp R - 

First-stage control mixing (longitudinal and vertical, lateral and directional) is 
followed by cumulative lateral stop limiting (of the authority of differential lateral 
and combined lateral inputs). Results of (vertical and lateral) second-stage mixing, 
©Fsp » ©ppp « ©RSP’ AND ©RPP are limited at the swashplate prior to driving the swivel- 
ing and pivoting upper-boost actuators. (In order that the swashplates move smoothly 
and not bind up, each is driven by a combination of swiveling and pivoting motions. 
Swashplate displacement is the sum of the two inputs.) The actuation dynamics are 
modeled as first-order lags, the outputs of which, ©jp R and Aj^p R * ma y inter- 
preted to be collective and lateral cyclic pitch angles represented in helicopter 
body axes, respectively. 

As described in reference 4, longitudinal cyclic pitch angle is scheduled with 
equivalent airspeed (fig. 34); actuation dynamics are modeled as a first-order lag. 
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Mechanical control 


system outputs are 

f A1CFRC , 
< BICFRC , 
l THOFRC , 


rotor hub collective 

AICRRC ] 

BICRRC V 
THORRC J 


and cyclic: positions 


which are passed to the ROTOR subroutine, 
in the CONTROL subroutine; table 8 gives 
cal flags. 


Table 7 is a summary of the variables used 
subroutine input/output variables and logi- 


Stability Augmentation System 


& w- is^sssi :rrxfi 


Ap 


static 

port 


= a.i) 




in (2y)sin(2S)q 


(57) 


where 


the angle between longitudinal axis and the static port line (52°) 
vnarrn r nro ccuva ( - fi K / 


q - the dynamic pressure (= (1 2 v r T 1 ^e static port line (52°) and 

cuUted to zer „ sld.sHp i6 ln 

38) 


angxe 18 S lven ln equation (58) where K. is a veloc 
dependent gain whose purpose is to wash out this rudder input at high R s P eeds (fig 


DRBYAW 6 


R. 


= (Ap in. H 7 0)(K 


equivalent 

pedal 


in. pedal 
in. H 2 0 


(58) 


Directional SAS yaw damping uses simple filtering with V = /n b «- u 

from a first-order lag in cascade with * i oa H-i g - Ve( l 40 knots ’ a change 

■— stssls 

d e rivauv:":a! ve n^ r “ !, 4S! 4 fzrs e ^?i"5 dyna " ic resp “ nse - 

4 m nZT^irt^ X'lZL, 


results ~or V 


eq 


- 7S and nn i ij-qj an a 40-4» give simil 

130 knots, respectively. More complete static and dynamic 
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model data may be found in volume II of this report, which gives the validation 
results. Included therein are static trim and stability derivative data as well as a 
summary of dynamic check results. 


Table 9 is a list of the SAS subroutine variables together with constants and 
conversion factors. Included is each variable, its FORTRAN mnemonic, its units, its 
common location if applicable, and its physical description. Table 10 is a list of 
the variables transferred between SAS and other subroutines. 


Electronic Control System 

Using the ECS of the CH— 47B, a researcher may either implement an experimental 
control system or, by designing explicit model-following laws, exercise the heli- 
copter's variable-stability capability. The ECS subroutine is a model of this system; 
subroutine inputs are the research pilot's cockpit control positions and the outputs 
are ECS signals sent to the mechanical controls subroutine, CONTROL. No specific ECS 
is documented in this report. It is anticipated that a particular ECS design will be 
developed along with an individual experiment, and will be documented at that time. 
However, during model validation, some simple procedures were developed which aid in 
properly linking the ECS to the rest of the model. A discussion of these follows. 

Prior to engaging the ECS (with flag IECSCON in subroutine CONTROL), the heli- 
copter is trimmed using the basic airframe and mechanical control system. In this 
case, the SAS must be turned off before trimming, since SAS inputs alter the cockpit 
control positions for trim. 


When the ECS is engaged the helicopter is flown by the research pilot; therefore, 
in the simulation the safety pilot’s inputs to the mechanical control system are 
disconnected as the ECS is turned on (see the CONTROL schematic, fig. 32). To avoid 
destroying the trimmed condition of the helicopter when the safety pilot's controls 
are disconnected, each trim cockpit control position is used as a bias which is added 
to the appropriate ECS input (which is zero at trim, by definition); this is shown in 
equations (59)-(62). 


where 


DLATTOT = DLATECS + DATOTIC 
DLONTOT = DLONECS + DBTOTIC 
DYAWTOT = DYAWECS + DRTOTIC 
DCOLTOT - DCOLECS + DCF’ TIC 


(59) 

(60) 
(61) 
(62) 


i 


r DAT0TIC > 


r DLATP > 

DBTOTIC 


DLONP 

> 

DRTOTIC 


DYAWP 

^DCTOTlCy 


J)C0LP^ 


trim 


Table 11 gives the ECS subroutine input/output variable definiti 


on 
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Slung Load 


Subroutine SLING models a baseline, externally suspended load in three degrees of 
freedom. This is accomplished by introducing three new state variables, each defined 
as a relative displacement of the load and helicopter body reference frames. Addi- 
tionally, terms which represent the effect of the slung load motion on the helicopter 
response are computed and passed to subroutine SMART. Simulation of the slung-load 
dynamics is optional and may be selected with flag ISLING. 

Figure 49 (taken from ref. 1) illustrates the geometry of the slung load, its 
attachment, and position relative to the helicopter. The baseline load data, which 
are included in the simulation model, is a "MIL-VAN" weighing 7500 lb. It is sus- 
pended on cables from tandem attachment points on the fuselage equally spaced about 
the nelicop' er c.g. It has been assumed that these attachment points may transmit no 
moments between the load and the helicopter. Referring to the figure: u L , X l , and 

V L are defined to be the longitudinal and lateral cable sway angles and the lateral 
differential cable angle, respectively. 

To compute slung-load aerodynamic quantities, velocities in the helicopter body 
reference frame at the slung-load c.g. are computed via equations (63)-(b3). 


USL 

U SL = U B + 

(l l 

+ V q B + l iA 

(63) 

VSL 

V SL = V B " 

(l l 

+ Vb - Vi. 

(64) 

WSL 

W SL = W B 



(65; 


Slung-load dynamic pressure, sideslip angle, and angle of attack, respectively, 
are computed in equations (66 ) — (68) . 


SQSL 

q SL 

= 1 P(U SL + V SL + W SL )1/2 

(66) 

BETSL 

^SL 

„ fSL\ 

= ' " L 

(67) 

ALFSL 

°SL 

/ w sl\ 

= arctan + 0 

W, / SL 

(68) 


Slung-load drag, sideforce, and yawing moment, respectively, are found from fig- 
ures 50-52 as a function of load angle of attack and sideslip angle. These data, 
normalized in the simulation model by load dynamic pressure, are taken from wind- 
tunnel tests. Prior to their use in the cable angle calculations, the load aerody- 
namic quantities ar* 5 resolved ir.tc '"he helicopter body reference frame, as in equa- 
tions (69)-(7 1) . 


XAEKSL 



-q 


SL 




(69) 
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(70) 

(71) 


*SL 


as inputs, suspension-cable angular accelerations 
omputed with the nonlinear second-order differential equations (72)-(74), from 


are c 
reference 1 


a AERl 

AMULrD ii T = “T 
L m L L L 




q B + l » 2 

\\ l l 


J L L L + V 




rp 


r v\ + ' Vl 1 (sin 0 + sin V " V 

Ij l 


(72) 


where 




[ (m L g)2 + X AERj] 
__ 


l/2 


-Y 


ALMLDD X T = 


aer l v b + 


l l + V 


L m L L L L I 


L. r B l_ T 2 L t | r B q B Vl 


( m L L L 


+ r B u L + (sin 4> + sin 

Ij Ij 


(73) 




ANULDD Vj = j 


"hj 8 ^ 

- " u~r cos 0 cos * V L " K v V L 

1 ; Ij 


(74) 


(During model validation, the value of K-, was changed from the original value of 
+1.8 to -0.03 to match BV dynamic-response data.) 

Integration results in cable angular velocities as in equations (75)— (77) . 


AMULD \ " J 

k 

dt 

(75) 

ALMLD \ = J 

k 

dt 

(76) 

ANULD \> T = J 

k 

dt 

(77, 
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and cable positions as in equations (78)— (80 ) : 


AMUL 

i ■ j 

) \ dt 

(78) 

ALML 

5 - J 

K dt 

(79) 

ANUL 


fi dt 

(80) 

At a straight and level flight condition 
living equations (72), (73), and (74) at 

, values of p^ , and may be 

steady state. Resulting trim values 

found 

are 


given in equations (81)-(83). 


AMULIC 

p = sin 1 

L I.C. 

-sin ( 

X AEBL 
D + 1 

(81) 

ALML I C 

X = -4> 

I.C. 



(82) 

ANULIC 

v L = 0 

I.C. 



(83) 


where I.C. represents the initial flight condition. By selecting flag ISLTRM, 
initial values of the slung-load states are computed at the same time that rhe heli- 
copter is being trimmed. 

In the original BV simulation model, the helicopter and slung load were modeled 
together as a coupled nine degree-of-f reedom system. However, since subroutine SMART 
was designed to handle only six degrees of freedom, the ARC model is somewhat modified 
from the original. Equations (84) -(89) are the nine degree-of-f reedom helicopter 
equations of motion in the helicopter body reference frame (ref. 1), where the under- 
lined terms are those which arise specifically from the slung load. 




• ‘AERO , n , , ,7 . x “L 

U " = — VB - 8 sin 0 + r B v B - jj- ( c i B w B - sin p l ) - ^ r B P B 


H 


4 

4 *h 


(84) 


^AFRO m L L . 

% = lf~ + 8 Sln * COS 0 + P B W B - r B U B + M" p B W B “ HT (r B q B + Vli 
H H L H 


% 


g sin X L 


( 85 ) 
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• M AERO ^XZ,, 2 / I 2Z " I xx \ R,J t 

q B " i~ + — <*S - pi> + Pn - ^ L 

yy yy \ x 


ORjc*rr ■ . 

OF POOH QU ALU'S 


* 
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Table 12 gives the SLING subroutine variable definition; table 13 is a list of 
subroutine input /output variables and logical flags. 


OPERATIONAL CONSIDERATIONS 


Real-time piloted simulation using a simulator cab and visual display requires 
the constant input information described in table 14. 

Additionally, in order that a pilot may land the helicopter model, a simple gear 
model has been devised. The landing gear subroutine is net actually executed; rather, 
subroutine BLAND has been modified so that the ground is contacted artificially (i.e., 
the gear reaction force is prescribed to be equal to the aircraft weight) and no 
reactive moments are calculated) . 


CONCLUSIONS 


A mathematical simulation model of the ARC CH-A7B helicopter has been purchased 
from the Boeing Vertol Company and implemented on the ARC Sigma IX computer. 

Volume I of this report includes engineering explanations of each model subroutine; 
also given are the appropriate assumptions and simplifications necessary to ensure 
the validity of a particular experiment. 


Volume II of this report gives a comparison among ARC and BV model dynamic 
response data and flight test data, together with ARC static-trim and stability- 
derivative data. Successful validation of the ARC model has been completed against 
BV model data. As with all mathematical models of physical systems, however, this 
model is not a perfect replication of the CH-47B helicopter. This is particularly 
true with a quasi-steady rotor dynamics model, the type implemented herein. To repre- 
sent specific aspects of the helicopter response more closely and to meet the needs 
of a particular simulation experiment, it may be desirable to modify the model 
described in this report. 
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APPENDIX A: FLAPPING AND CONING EQUATIONS 


Using Wheatley-Bailey theory (refs. 6 and 7), flapping and coning angles are 
computed in the shaft-normal-plane-wind reference frame. Due to pitch-flap coupling 
(6 3 ), the solution for coning and flapping angles (a^a^bi) is coupled with the 
definition of swashplate cyclic and collective pitch angles (A x .Bj ,O 0 ) , as shown in 
equations (A1)-(A7). Additionally, coning angle is a function of rotor thrust, 
defined in equation (Al). 
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AORR 
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AICFR 
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BICRR 
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(A4) 

(A5) 

(A6) 

(A7 ) 


The purpose of this appendix is to provide the algebraic steps necessary to 
decouple these equations, eventually resulting in the model equation (18). i\,llowi:v, 
is the step by step decoupling of the equations, reproduced from reference 4. 


1. Substituting for (2C T /ao) in the a„ equation: 
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2. Substitution for O 0 and B x : 
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. K 

jacRg S 


(I * I ') 


+ a j (0 . 0) + b 1 (2 V iK ) 


(All) 


= u'^l + | ,) 


2yBl + 2 A + 0 


“2 


<*0 + u ‘) 


(A12) 


3. After defining coefficients as indicated above, the equation has the form: 

Aa 0 4- Ba 1 4- Cb 1 = J (A13) 

4. Rearranging the a 1 equation: 


&-£)•>- 


n* + — 

2 3 


2v-0„ yC fc „ 3 u 2 B^ B^ 16q F>R [l + ( U 2 / 2 ) ] [ (1/4) - (p 2 /8) ] 


0 + " tv 


8 


(pacRgJO /1 F>R 


5. Substituting for O and B, 


(A 14 ) 


(i " l) a i 


(o; + Vo) 3 


+ ^ _ 

2 2 


( B ^c ? + v.)(i + ¥) 

U °tw 16 q F,R [1 + (u 2 / 2 ) ] t Cl/ 4 ) - (y 2 / 8 )] 


(- f v) + (1 - t) + >’f e (i + ¥) 


(pacRjn)/I p>R 


D 
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t!± _ B . (i + 3 y 2 \ 

3 Bl °2 v V T 


F 

yo 


tw Ife q FiR [i + (y 2 /2) 1 [ (1/4) - (y 2 /8) i 


(pacRgft) /Ip^ R 


(A 15 ) 


(A 16 ) 
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’ "* p,J 


7T 




6. The definition of the above coefficients results in the form of equa- 
tion (A17): 


Da 0 + Ea 1 + Fb x = K 


7. Substituting for A lc in the b x equation: 


b, = 


4ua, 


3 [ 1 + (u 2 /2) ] 


+ A,' + K a 


e i 


16p F,R [ 1 " (M 2 /2)] 

(p.cR*n)/I 


CA 17) 


(A18) 


l- 
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-4 u ... 16 Pv vt 1 " C U 2 /2 ) ] 

*0 : + a : (-K ) + bjd.O) = A[ ^ 

3[1 + (p /2) ] c 2 (pacRgfi) /Ip R 


G 


(A19) 


H I L 

8. After the above definitions, the equation has the form: 

Ga 0 + Ha x + Ib x = L (A20) 

ln the equations (A13) , (A17), and (A20) , which are the same as the 

text matrix equation (18), are solved for a 0 , a x> and b x . 
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APPENDIX B: INFLOW DYNAMICS SOLUTION 


In the original version of the model (developed by Boeing Vertol Company) from 
which this model was adapted, the inflow ratio was modeled by the equation (Bl) 
expression, including a first-order lag (ref. 1). Past cycle values of C Tf>r , Af,R, 

and y F R were used, so no iteration on the current value of y F>R was performed 
using this implementation. 


ALAMFR , 
ALAMRR A F,R 


\ ' 

A F,R 


Ct f,r , F RF(FR) t r,f 

1 

_ 2(p F,R + X F,R } 1 2(P R,F + X R,F ) _ 

• s + 1 

A f,r J 


(Bl) 


where 


w 


F,R 


w 


F,R 


F,R " R B (fi’ iR - r F>R ) R B ^ ^ f ,R 


u f,r 


F,R 


as defined in model equation (12). 

A more exact real-time solution was obtained by Boris Voh , who represented the 
above as a differential equation and solved it using a local linearization method 
implemented as subroutine LOLIN (ref. 12). Following is the solution method, using 
the forward rotor equation as an example: 


X F = 


W F 

C T D F C T 

l T RF R 
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R B F F 

2(P F + A F ) l/2 2(y| + 

r . s + 1 
_ X F 


{v + 1 )’A(v + 1 ) ■ 
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RF R 


2(ya, + Aj) l/2 2 (p R + A r ) 1/2 


(B2) 
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W T 
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F R Q r 

BpF \ 
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F 1 " F > 
C r 


2 K + a r> 


l/2 


°F C T 
RF R 


V“f 2(y* + A*) l/2 


2(4 + 4> 1/: 


(B4) 


(B5) 


“ F * “ vM F ^ 'F / " vm R ' R' 

Following are the definitions necessary for the application of LOLIN to this problem: 
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Description 


LOLIN 

definition 


Fngineering 

definition 


Nonlinear function 


FN 



Partial derivative of nonlinear FT 

function with respect to time 



Jacobian of system 


3A f 

3A t , 

r 

3A f 

3X r 
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3A I; 
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3A t 
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System state vector 


SI 



where 



is defined above in equation (B5) , the time derivative of the function, 


tm « 


r i 

K 


0 

F 

_ 

5L 


0 
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m mt 




and the four elements of the Jacobian may be calculated as in equations (B6)-(B9) 
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3X F " ' T A, 


a F°F 


c t f x f 
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(1 - sin 3 


RF 
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)Ad* 


3A. 


|-(A r / Ur ) - o'. 25 ~ 


RF 


(B6) 
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where 
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3a f x A. 
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T F *FR F 
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FR 
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FR A F 
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TABLE 1A.- CONCLUDED. 



* CH(77) Rear rotor hub vertical force, transformed to helicopter 

body reference frame 



TABLE IB.- ROTOR CONSTANTS AND CONVERSION FACTORS 




.1042 Forward rotor hub mean aerodynamic chord 


Common M . .. 

Simulation Engineering Units location (if ora i naJ - Physical description 

mnemonic variable applicable) Value 



TABLE 2.- ROTOR SUBROUTINE VARIABLE DEFINITION. 


Input variables 


Output variables 


Variable 

Common 

location 

Subroutine 
of origin 

Variable 

Common 

location 

Subroutine of 
destination 

AICFRC 

CH(39) 

CONTROL 

ALARFR 

CH(52) 

AERO 

A[CRRC 

CH(38) 

1 

ALARRR 

CH(53) 



BICFRC 

CH(37) 

1 

AMARFR 

CH(54) 



BICRRC 

GH(36) 

▼ 

AMARRR 

CH(55) 



HCG 

A(176) 

SMART 

ANARFR 

CH(56) 



OMEGPF 

CH(115) 

ENGINE 

ANARRR 

CH(57) 

1 


OMEGPR 

CH(116) 

ENGINE 

QAERFR 

CH(64) 

ENGINE 

PB 

A(37) 

SMART 

QAERRR 

CH(65) 

ENGINE 

QB 

A(38) 

SMART 

TMGN01 

CH(100) 

AERO 

QGOVFR 

CH(257) 

ENGINE 

TMGN05 

CH(104) 



QGOVRR 

CH(258) 

ENGINE 

TMGN06 

CH(105) 



RB 

A(39) 

SMART 

XAERFR 

CH(72) 



SBETFS 

CH(50) 

AERO 

XAERRR 

CH(75) 



THOFRC 

CH(42) 

CONTROL 

YAERFR 

CH(73) 



THORRC 

CH(43) 

CONTROL 

YAERRR 

CH(76) 



UB 

A(58) 

SMART 

ZAERFR 

CH(74) 



VB 

A(39) 

SMART 

ZAERRR 

CH(77) 


r 

WB 

A(60) 

SMART 







Lo 

1 « 

gical flags 


Required 

input data 

Flag 

Common 

location 

Function 

Variable 

Common 

location 

Description 

NGREFF 

ICH(7) 

Ground effect correction 

DXCG 

CH(68) 

Position of actual 



of thrust off/on (0/1) 

DYCG 

CH(69) 

helicopter c.g. 

N STALL 

ICH(5) 

Rotor stall modification 

DZCG 

CH(70) 

relative to its 



of thrust and torque 



reference (fig. 30) 



off/on (0/1) 




NTRQCR 

ICH(6) 

Empirical correction 






of rotor torque 






off/on (0/1) 






TABLE 3A.- AERO SUBROUTINE VARIABLE DEFINITION 
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TABLE 3A. - CONTINUED. 




T EMA 1.689 p/p CH(92) Conversion factor between TAS (kt) and CAS (ft/sec) 
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TABLE 3A. - CONCLUDED. 
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TABLE 4 . - AERO SUBROUTINE TRANSFER VARIABLES, INPUT DATA AND 

LOGICAL FLAGS 


Input variables 

Output variables 

Variable 

Common 

location 

Subroutine 
of origin 

Variable 

Common 

location 

Subroutine of 
destination 

ALARFR 

CH(52) 

ROTOR 

BETAFS 

CH(59) 

SAS 

ALARRR 

CH(53) 



FAX 

A(136) 

SMART 

AMARFR 

CH( 54) 



FAY 

A(137) 



AMARRR 

CH(55) 



FAZ 

A(138) 



ANARFR 

CH(56) 



TAL 

A(155) 



ANARRR 

CH(57) 

1 

f 

TAM 

A(156) 


r 

TEMA 

CH(92) 

SMART 

TAN 

A(157) 



TMGNOl 

CH(100) 

ROTOR 




TMGN05 

CH(104) 

ROTOR 





UB 

A(58) 

SMART 





VB 

A(59) 

SMART 





WB 

A(60) 

SMART 





XAERFR 

CH(72) 

ROTOR 





XAERRR 

CH(75) 







YAERFR 

CH(73) 







YAERRR 

CH ( 7 6 ) 







ZAERFR 

CH( 74) 







ZAERRR 

CH (77) 
n 

1 

r 



i 



Required Input Data 


Variable 

Common 

location 

Description 

DXCG 

DYCG 

DZCG 

CH(68) 

CH(69) 

CH(70) 

Position of actual 
helicopter c.g. relative 
to its reference (fig. 30). 

L - ■ - - 


60 



TABLE 5A.- ENGINE SUBROUTINE VARIABLES 



TABLE 5A. - CONTINUED. 









TABLE 5B.- ENGINE CONSTANTS AND CONVERSION FACTORS 








Logical flags 


Flag 

Common 

location 

Function 

I STEADY 

ICH(4) 

Zeros n after rigid body 
states have been trimmed 
off/on (0/1) 


TABLE 6.- ENGINE SUBROUTINE TRANSFER VARIABLES 



Input variables 

i 

, Output variables 

Variable 

Common 

location 

Subroutine 
of origin 

| Variable 

Common 

location 

Subroutine of 
destination 

DCOLTOT 

CH(210) 

CONTROL 

| OMEGPF 

CH(115) 

ROTOR 




, OMEGPR 

CH(11 6) 

ROTOR 

IBEEPl 


Simulator cab 




IBEEP12 


Simulator cab 

| QGOVFR 

CH(257) 

ROTOR 




, QGOVRR 

CH(258) 

ROTOR 

QAERFR 

CH(64) 

ROTOR 




QAERRR 

CH(65) 

ROTOR 

1 

1 






TABLE 7A. - CONTROL SUBROUTINE VARIABLE DEFINITION 


0 



DAICRRC A, deg Rear rotor lateral cyclic pitch, body reference frame 






TABLE 7A. - CONTINUED. 
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TABLE 8.- CONTROL SUBROUTINE TRANSFER VARIABLES 


Input variables ' 

Variable 

Common 

location 

Subroutine of 1 
origin 1 

DCOLECS 

CH(275) 

ECS ' 

DCOLP 

CH(206) 

Simulator cab 

DLATECS 

CH(272) 

ECS j 

DLATP 

CH(203) 

Simulator cab , 

DLATSAS 

CH(17) 

SAS , 

DLONECS 

CH(273) 

ECS 

DLOMP 

CH(204) 

Simulator cab . 

DLONSAS 

CH(18) 

SAS , 

DYAWECS 

CH(274) 

ECS , 

DYAWP 

CH( 205) 

Simulator cab , 

DYAWSAS 

CH(19) 

SAS , 

IAND 

IA(29) 

Simulator cab , 

IANU 

IA(30) 

Simulator cab , 

ILWD 

IA(33) 

Simulator cab , 

IRWD 

IA(34) 

Simulator cab , 

VEQ 

1 

A(75) 

SMART , 


Output variables 


Variable 


AICFRC 

AICRRC 

B1CFRC 

BICRRC 

THOFRC 

THORRC 


Common 

location 


CH(39) 

CH(38) 

CH(37) 

CH(36) 

CH(42) 

CH(43) 


Subroutine of 
destination 


ROTOR 

ROTOR 

ROTOR 

ROTOR 

ROTOR 

ROTOR 


Logical flags 


Flag 

Common 

location 

Function 

IDCPT 

ICH(3) 

Differential collective pitch 
trim off/on (0/1) 

IECSCON 

ICH(2) 

Electronic control system off/on (0/1) 

IMHIS 

— 

Simulator cab off/on (0/1) 

RSASP 

CH(282) 

Lateral SAS off/on (0/1) 

RSASQ 

CH(283) 

Longitudinal SAS off/on (0/1) 

RSASR 

CH(284) 

Directional SAS off/on (0/1) 
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TABLE 9B. - SAS CONSTANTS AND CONVERSION FACTORS 
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TABLE 11.- ECS SUBROUTINE TRANSFER VARIABLES. 


r — — 

Input variables i Output variables 

Variable 

Common 

location 

Subroutine of „ . , _ 

. . i Variable 

origin 

Common 

location 

Subroutine of 
destination 

DCOLP 

DLATP 

DLONP 

DYAWP 

CH(106) 
CH(203) 
CH ( 204 ) 
CH(205) 

1 ' f " 

Simulator cab ■ DCOLECS 
Simulator cab i DLATECS 
Simulator cab • DLONECS 
Simulator cab ' DYAWECS 

CH ( 2 7 5 ) 
CH(272) 
CH(273) 
CK(274) 

CONTROL 

CONTROL 

CONTROL 

CONTROL 
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TABLE 12B. - SLING CONSTANTS AND CONVERSION FACTORS 
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TABLE 13.- SLING SUBROUTINE TRANSFER VARIABLES, INPUT DATA AND 

LOGICAL FLAGS. 


Input variables i Output variables 


Variabl e 

Common 

location 

Subroutine^ Varlable 
ot origin 

Common 

location 

Subroutine of 
destination 

CPHI 

A( 11) 

SMART ' PBDS 1 

CH( 251) 

SMART 

CTHT 

A(13) 

1 QBDS 

CH(252) 


PB 

A(37) 

1 RBDS 

CH(253) 


PBD 

A(35) 

1 UBDS 

CH(248) 


PHIR 

A(4) 

1 VBDS 

CH(249) 


QB 

A(38) 

1 WBDS 

CH(250) 


QBD 

A(56) 

1 



RB 

A(39) 




RBD 

A(57) 




RH02 

CH(101) 




SPHI 

A(10) 

1 



STHT 

A(12) 

1 

« 



UB 

A(58) 

1 

| 



UBD 

A(413) 

1 

i 


VB 

A(59) 

1 

1 



VBD 

A(414) 

1 

| 


l 

WB 

A(60) 

1 



XIXX 

A(116) 

1 

j 



XIXZ 

A(119) 

| 



XIYY 

A(117) 

J 



XIZZ 

A(118) 

j 



XMASS 

A(130) 


.. 

J p-J 


Required Input Data 


Variable 

Common 

location 

Description 

BJSL 

CH(266) 

Moment of inertia of slung load about 
load vertical axis 

BLSL 

CH(267) 

Average cable length below attachment 
point 

BRSL 

CH(2A8) 

Vertical distance between hook attach' 
ment point and aircraft c.g. 

SASL 

SMSLIC 

CH(297) 

Cable separation distance 
Slung load mass 

THESL 

CH(269) 

Angle between load x-axis and heli- 
copter x-axis 

WGHTSL 


Slung load weight 


Logical Flags 


Flag 

Common location 

^ Function 

ISLING 

ISLTRM 

1CH(1) 1 

ICH(8) 

Slung load subroutine option off/on (0/1) 

Slung load trim in straight level flight off/on 
(0/1) 



TABLE 14.- REQUIRED INPUT DATA FOR OPERATIONAL SIMULATIONS 


Variable 

Common 

location 

Units 

Physical description 

DXCG 

CH(68) 

in . 

Position of actual helicopter 

DYCG 

CH(69) 

in. 

c.g. relative to its refer- 

DZCG 

CH(70) 

in * 

ence (fig. 30) 

WAITIC 

A (24 2) 

!b f 

Helicopter weight 

XIXXIC 

A(243) 

slug-ft 2 

Helicopter moments and prod- 

XIYYIC 

A(244 ) 

slug-ft 2 

uct of inertia 

XIZZIC 

A(245) 

slug-ft 2 


XIXZIC 

A(246) 

slug-ft 2 


XP 

A(171) 

ft 

Position of pilot, in heli- 

YP 

A(172) 

ft 

copter body axes, relative 

ZP 

A(173) 

ft 

to c.g. of aircraft 


i 







ORIGINAL PAGE IS 
OF POOR QUALITY 
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ABRO | 

J 


FROM KNQINC 
ENGINE [ OOOVFR i 
TORQUE | QGOVRR j 


(FROM 5) 


^AVAILABLE j 


ISBITF I 

icr* .'ll 


SLFR 

SOFH 

SHF* 


(FROM 1)- 


SLRR 

SDRR 

SHRR 


17 COMPUTE ROTOR 
AERODYNAMIC 
MOMENTS DUE 
TO HUB FORCES 


ALFR1 

AMFR1 

ANFR1 


ALRR1 

AMRR1 

ANRR1 


XAERFR 
Y AERFR 
Z AERFR 


ALFR1 

AMFR1 

ANFR1 


ALRR1 

AMRR1 

AWRR1 


19 COMPUTE TOTAL 
ROTOR MOMENT S 
(HELICOPTER BOOY 
REFERENCE FRAME) 


ALARFR 

ALARriR 


AMARFR 

AMARRR 


ANAAFR 




ALFR2 

AMFR2 

ANFR2 


ALRR2 

AMRR2 

ANRR2 


I 


18 TRANSFORM ROTOR 
HUB MOMENTS TO 
HELICOPTER BODY 
REFERENCE FRAME 


ALFR2 

AMFR2 

ANFR2 


A L RR2 

AMRR2 

ANRR2 


XAKRRR 
I YAIRRR 
ZAERRR 


(StNIFRI 

ICOSIFRj 


ISIMIHRI 

ICOSIRRj 


ISISKTRRI 
jCBITRR | 


I TO 


1 

aero! 


16 TRANSFORM ROTOR 
AERODYNAMIC FORCES 
TO HELICOPTER BODY 
REFERENCE FRAME 
(THROUGH SHAFT 
INCIDENCE AND ROTOR 
SIDESLIP ANGLES). 


XAERFR 
Y AERFR 
Z AERFR 


XAERRR 
YAERRR 
ZAERRR | 


("from 

] SMART 

1 HELICOPTER 


LONGITUDINAL j 

BODY 1 


UB 


VELOCIT 
| REFERENlE I 
jJFRAME) 


3 COMPUTE 
ROTOR/ROTOR 

ALAMFR 

15. COMPUTE ROTOR 
INFLOW RATIO 


INTERFERENCE 

PARAMETER 

(REVERSE 

^ALAMRR 

ALAMFR , ALAMRR 

_ CTFR 

FUNCTIONS FOR 
REARWARD 

BOFRF 

(SUBROUTINE 
LOLIM CALLED) 

CTRR 

FLIGHT}. 

BOFRF , BDFFR 

BDFFR 

(FROM 13) 






5 


sin ip, cot «p, tin ip, cot tp 


IS i NIFR j 
tCOSIFR! 


sm.cot ip 


HFR 

YFR 

TFR 


HRR 

YRR 

TRR 


ISIHIFRi 

1COSIFRI 

sin, cot ip 


1 14. COMPUTE 

n 

BMWFR 

ROTOR BLADE 

A. ROTOR HUB 


BMWRR 

MASS MOMENT 

PITCHING MOMENT 


AMHBFR , AMHBRR 

_ SBFR 

NUMBER OF 

B, ROTOR HUB 


SBRR 

BLADES PER HUB 

ROLLING MOMENT 

ALHBFR , ALH8RR 

_ SEFR 

FLAPPING HINGE 

SERR 

OFFSET 

iAIFR! 

A1 RR 


[BIFRl 

[BIRR] 

I 


FROM 

SMART 

HELICOPTER 

ANGULAR 

VELOCITIES 

(BOOY 

REFERENCE 

FRAME) 


SINES AND 
COSINES 
OF SHAFT 
INCIDENCE 
ANGLES, 

*F. *P 


TRANSFORM AIRFRAME 
ANGULAR RATES TO 
SHAFT NORMAL PLANE 
WIND REFERENCE FRAME 
(THROUGH SHAFT 
INCIDENCE AND ROTOR 
SIDESLIP ANGLES ) 


PFR] 

| PRR 

QFR 

QRR 

rfr! 

, [RRR 


ISBCTFRI 

I cbetfr! 


(SBETRR 

ICBETRRf 


(FROM S) 


I FROM: "1 

| ENGINE | 

u:^_j 

loME GPF 

OMEGPR 


ROTOR 

BLADE 


RBFR 


RADIUS RBRR 


AIR RHO 
DENSITY 


RFR 


RRR 


(FROM 12) 


TAN ( -. 3 ) 


COMPUTE 

A. ROTOR ANGULAR 
VELOCITY CORRECTED 
FOR HELICOPTER 
YAW RATE. 

OMEGFR , OMEGRR 


KBETF KBETR 


b ROTOR TIP SPEED 

VTIPFR , VTIPRR 


C. ROTOR ADVANCE 
RATIO 

AMUFR , AMURR 


D INFLOW RATIO 

(FREESTREAM 

COMPONENT) 

ALMPFR , ALMPRR 


j U FR 
IWFR 


. (FROM 4) 


URR 

WRR 


jAICRR 

Ibicrr 

[thorr 


e. COMF 


11. CORRECT 

CONTROL INPUTS 
FOR PITCH-FLAP 
COUPLING 


AICRF 


AICRR 

BICFR 


b ICRR 

THOFR 

, 

THORR 


TAICFR2 

BICFR2 


IS I NIFR! 
ICOSIFRj 


ISINIRR) 

COSIRRj 


COMPUTE VELOCITY 
IN ROTOR SHAFT 
NORMAL PLANE 
(RESOLVED 
THROUGH SHAFT 
INCIDENCE ANGLE) 


UFR2 
V FR2 
WFR2 


URR2 
V RR2 
,[WRR2 


r ” 


FROM: 

! SMART 

I HELICOPTER 
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igure 3.- Helicopter rotor center of gravity positions relative to rotorcraft 

center of gravity (ref. 1). 


Figure 4 


^actual 



- Actual versus baseline helicopter center of gravity position. 
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Figure 5.- Reference frame transformation through shaft incidence angles 



Figure 6.- Reference frame transformation through rotor sideslip angles. 
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Figure 9.- Rotor stall thrust 


coefficient correction (subroutine RSTALL) 
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Figure 12.- Empirical correction of rotor torque coefficient 

(in-line calculation) . 



Figure 13. Rotor-on-rotor interference terms: d^, (forward flight) and 

and d' (rearward flight) 
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Figure 15.- Rotor on rotor interference term. 
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Figure 16.- Fuselage drag data (table: FDOQT) . 
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Figure 20. ~ Fuselage pitching moment data (table: FMOQT) 
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Figure 21.- Fuselage yawing moment data (table: FNOQT) . 



Figure 22.- Actual helicopter versus wind 
tunnel model center of gravity. 
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Figure 23.- Engine, governor, and shaft dynamic s block diagram. 
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Figure 26.- Power dynamics time constant (table: POWT, TAUPT) 



Figure 27.- Variable power limits (table: POWERT, E27T). 
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Figure 30.- Percent topping power (table: OMEGAT, POMEGT) 
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Figure 33.- Longitudinal cyclic position gradient stabilization with 
differential collective pitch trim. 
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Figure 35.- Longitudinal stability augmentation system. 
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Figure 36.- Lateral stability augmentation system. 
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Figure 37.- Directional stability augmentation system 
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Figure 44.- Lateral axis dynamic response SAS OFF and ON; 
= 75 knots, weight = 33,000 lb, nominal c.g. position. 
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Figure 50.- Slung load drag force (table: 3LDQT). 
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ure 51.- Slung load side force (table: SLYQT) . 
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Figure 52.- Slung-load yawing moment (table: SI.NQT) . 



